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by 

A. H. Pacella** and F. Erdogan 
Lehigh University, Bethlehem, Pa. 


Abstract . Using the filament model developed in the previous 
paper, the elastostatic interaction problem between a penny-shaped 
crack and a slender inclusion or filament in an el as ti c ma tr i x is 
formulated. For a single filament as well as multiple identical 
filaments located symmet ri cal 1 y around the crack the problem is 
shown to reduce to a singular integral equation. The solution of 
the problem is obtained for various geometries and filament-to- 
matrix stiffness ratios, and the results relating to the angular 
variation of the stress intensity factor and the maximum filament 
stress are presented. 


1 . INTRODUCTION 

In studies relating to the fracture of filament-reinforced 
composites and reinforced concrete in the presence of internal 
cracks, in practice two of the most important factors one needs 
to evaluate are the effect of the reinforcements on the stress 
state around the leading edge of the crack and the effect of the 
crack on the maximum stresses in the reinforcements. Basically 
the problem is that of interaction between an internal crack in 
an elastic matrix and elastic inclusions of different thermo- 
elastic properties and of given size, shape, orientation, and 

* This work was supported by The National Science Foundation 
under the Grant GK-11977 and by the National Aeronautics and 
Space Administration under the Grant NGR-39-007-011. 

**Present Address: Bettis Atomic Power Lab., Westinghouse 
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distribution which are imbedded in the matrix. The general prob- 
lem is clearly intractable. However, if the medium fulfills 
certain conditions of symmetry with respect to the geometry and 
the external loads, under some simplifying assumptions it is 
possible to solve the problem within a reasonable degree of 
accuracy and with a reasonable amount of computational effort. 
Such a solution will be given in this paper. 

Specifically, the following are the symmetry conditions 
which are assumed to be satisfied: 

(a) the crack is a plane, circular (i.e., penny-shaped) 
internal crack; 

(b) the filaments are circular cylinders imbedded in the 
matrix with their axes perpendicular to the plane of the crack, 
z = 0 ( Fi gure 1 ) ; 

(c) z~0 is a plane of symmetry with respect to the geom- 
etry of the composite medium and the external loads; and 

(d) the overall dimensions of the matrix are large in com- 
parison with the local dimensions of the crack-filament region 
(i.e., the crack diameter, 2a, the distance of the ith filament 
from the crack center, b^- , and the filament length, 2c.) so that 
in the analysis the matrix may be treated as being infinite. 

The primary assumption to render the problem analytically 
tractable is that the Poisson's ratios of the filaments be equal 
to that of the matrix. Under this assumption it was shown in [1] 
that if the Green's functions (for a concentrated body force) of 
the matrix are available, the problem may be reduced to a system 
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of (two-dimensional) integral equations for layers of (fictitious) 
body forces imbedded along the filament- matrix interfaces. The 
following simplifying assumptions are then made in this paper in 
order to keep the analysis and the numerical work to solve the 
problem within manageable proportions: 

(i) The filaments are identical, i.e., c-j = ... = Cj^=c, 

ri= ... = r^=ro, = . . . = E„ = . v, = . . . = = v. where 2c 

is the length, r^ is the radius, is the Young's modulus, is 
the Poisson's ratio, and N is the number of filaments, and v is 
the Poisson's ratio of the matrix (Figure 1); 

(ii) c/r^ is sufficiently large so that a one-dimensional 
treatment of the auxiliary inclusion (having the elastic constants 
E^- E and v, E being the modulus of the matrix) is justified (see 
[ 1 ]); 

(iii) The filaments are distributed symmetrically (i.e., 
with equal angular spacings) on a circle in z-O plane concentric 
with the penny-shaped crack (i.e., b-| = ... =b|^ = b), they are 
parallel to the z axis, and their mid-points lie on z=0 plane; 

(iv) The filament radius r^ is sufficiently small in com- 
parison with other dimensions of the medium so that the effect of 
the radial component R of the (fictitious) body force acting on 
the filament-matrix interface and the effect of the angular varia- 
tion of the axial component Z of the body force may be neglected. 
It was shown in [1] that as a result of this assumption the effect 
of the. end tractions (or the body forces distributed at the ends 
z=+c) may also be neglected; and 

(v) The external load is the uniaxial stress cr^z ~ ^o 
to the matrix away from the crack-filament region. 
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From the viewpoint of the elasticity solution, clearly the 
most drastic assumption is that regarding the distribution of the 
body forces. However, the results and some comparisons given in 
[1] indicate that the error introduced into the solution as a 
result of this assumption would be within an acceptable range. 

The physical problem described in this paper was also consid- 
ered in [2], where the problem was assumed to be essentially 
axisymmetric and the filaments were replaced by a discrete set of 
concentrated body forces acting in z-direction. 

2. FORMULATION OF THE PROBLEM 

Consider now the problem described in Figure 1 subject to 
the symmetry conditions (a)-(d) and the simplifying assumptions 
(i)-(v) outlined in the previous Section. Since the filaments 
are identical and symmetr i cal ly located, the body forces Z(z) 
imbedded in different filament-matrix interfaces will be the same. 
Hence if one can obtain the Green's functions for a cracked 
matrix, it is possible to formulate the problem in terms of a 
single integral equation with Z(z) as the unknown function. This 
integral equation may be obtained by matching either the z-compo- 
nents of the displacement vector or the strains in the cracked 

matrix subjected to the external loads and Z(z) and that in 
the auxiliary filaments subjected to the external loads -Z{z). 

In using the equality of the strains in the derivation the single- 
valuedness condition of the displacements must separately be 
satisfied. The solution for the cracked matrix may be obtained 
by superimposing four solutions described in Figure 2. Note that 
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there are no cracks In problems described by Figures 2a and 2c. 

Let the filament 1 be located on the x axis (i.e., at r=b, 

e = 0) and let (6,<f),z) be the cylindrical coordinates referred to 

the center of filament 1. Let ^ 2 z’*’’^zz strains corre- 

f 

spending to problems described by Figures 2a,..,2d, and be the 

axial strain in the auxiliary filament. From [1] it follows that 


£ +e +e +£ 

zz zz zz zz zz' 


( 5 = r^ , (j) = , -c < z < c) , (1 ) 

where r^ is the radius, 2c is the length of the filament, and 4 >q 
is an appropriate angle describing the generator of the cylinder 
along which the strains are matched (Note that the body force 
Z(r^,(}),z) is assumed to be independent of 4>). The filament strain 
is given by [1] 


= - r„(E^-ET I 


( - c < z < c ) . 


( 2 ) 


Within the confines of the simplifying assumptions made in the 
previous section regarding the body forces R and Z, the strains 
and were obtained in [1] as follows: 


zz 


zz 


zz 


(3) 


zz 


(z) = + / ni22(z.t)Z(t)dt] , 

(-C < z < c ) , 


(4) 
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, ^'"o ,3K(k) 1 > ■ pQ ^(|,x 


H- !l°l 3E(k) ^ 


PQ 


E(k) 


N 


1 


+ irr^ ( t-z ) y — 9 0/0 

° i=2 [d^ + {t-z)2]3/2 


(1 - 2y 


, 3y(t-z)^ ' 

P o / > 

+ {t-z)^ 


4tt ( 1 - V ) 


1 


1 (1+v) (3 - 4v) ’ 


2 /, 2 ^ ,^2 
pQ "’Pq + ’ 


" 3 - 4v ’ 

k = — 9- . 


d. = b[2(l - cos , (i = 2.. . ,N) , (5) 

where K(k) and E(k) are complete elliptic integrals, N is the 
number of filaments, and b is the distance from the center of the 
crack to the midpoint of the filament {Figure 1). 


To obtain and the solutions of the crack problems 
shown in Figures 2b and2d (where the only external loads are the 
self-equilibrating crack surface tractions) are needed. In Figure 
2b the crack surface tractions are 

°zz = - Pq ’ Prz " ° = P0Z • 

This is an axisymmetric problem the solution of which is known 
(see, e.g., [4]). In the problem described by Figure 2d the crack 
surface tractions are equal in magnitude and opposite in sign to 
the stresses at z=0 plane obtained from the solution of the prob- 
lem given by Figure 2c. In Figure 2c, because of symmetry, the 
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shear stresses at z=0 plane will again be zero. Hence for the 
problem 2d the crack surface tractions will be 

= - p(r,9) , - 0 = . 

where p(r,6) is the normal stress at z=0 plane in Figure 2c. 


p(r,9) may easily be calculated in terms of Z(z) by using 

the basic solution given in [3] (see also [1]). However, since 

r <<c and r <<a (where r is the radius, 2c is the length of 
0 0 0 

the filament and a is the crack radius), to calculate p(r,9) for 
simplicity it will be assumed that the concentrated body forces 
2TirQZ(z) act along the lines (r = b, --c<z<c, 6^ = 27ri/N), 

i=l,2,,.,N (i.e., along the axes of the filaments), rather than 
Z(z) acting along the cylindrical surfaces ('5 = r^, - c < z < c , 

0_<(l)<2Tr). Thus, following [3] we find 

r c 

= 4/1^,^ / ^ -L (1 - 2v + ^-) dt . 

' ■'-c i = lp. p, 

1 1 

p? = - 2br cos(e--^) + . (8) 

For the solution of the crack problem it will be necessary to 
express p(r,0) in product form. Since p is periodic (and even) 
in 9 this may easily be done, giving 
“ « Nn 

p(r,0) = I b (^) cos(Nne) , (0<r<a, O£0<27r), (9) 

n = 0 " ® 


where 

1 a 2 tt 

b = — 2 / / p(r,0)rdrde , 
Tia^ 0 0 


^/M,T^a27r ^nN + 1 

bn "" ^ ^ cos(nN0)drd9 , 

0 0 

(n = 1,2,...). 


( 10 ) 
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Assuming that crack radius a is the length unit for the medium, 
without any loss in generality for the remainder of this paper 
we wi 1 1 1 et a ~ 1 . 


The crack problem described in Figure 2d may thus be consid- 
ered as a mixed boundary value problem in elasticity for a half 
space z>0 subject to the following boundary conditions: 

a^^(r,e,0) = - p(r,e) , (0<r<l, 0<e<2-ir), 

= 0 » (l<r<«’, O<0<2 tt), 

= 0 = aQ^(r,0,O) , (0<r<<», O<0<2 tt). 

(ll.a-c) 

Similar mixed boundary value problems have been treated in [5] - 
[7]. As shown in, for example, [8], because of (11. c) the prob- 
lem may be formulated in terms of a single harmonic function (j), 
which may be expressed as 


oo oo _ 

(|.(r,0,z) = I cos(nN0) / ^ f^(a)e“^^J^j^(ra)da , (1 

n = 0 0 

where f^ is the unknown function. Using (12) and the relations 


= 2(l-v) If - z 0 , 

2 3 

9 4 ) , 9 4 ) 

^zz ~ 2 " ^ 3 ’ 

az^ az'^ 


from (9), (11. a) and (11. b) it is found that 


/ af^(a) J^j^(ra)da = - ’ (0<r<l), 

0 

00 

/ = 0 > (l<r<«‘), (n = 0,l,2,..,). 

0 


(14) 
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After some manipulations the solution of the dual integral equa- 
tions (14) may be expressed as [4» 9] 




= 2 


-2nN 


nN /2nN + l 

M , 

k=0 \ k , 




n-k 


2n - 2k + 1 


(15) 


Substituting now (12) and (15) into (13. a) the desired strain 
may now be expressed as 

00 

e^^(r,e,z) = - ^ a^b^cos(nN0)[(l - 2v) ( r , | z | ) 

+ |z|G^(r,lz|)] , (|z| >0), (16) 

where 

00 1/0 

F^(r,z) = / a ’ 

00 0/0 

G^(r,z) = / ot^^^Jj^l^(ra) ’ (z>0)« (17. a, b) 

0 

The functions and G^ may be evaluated by using Filon's inte- 
gration formula [9]. However, a simpler way to evaluate these 
integrals would be to use the series expansions given in [10, 
p. 400]*. 

The problem described by Figure 2b is a special case of 
Figure 2d, replacing p(r,e) in (11. a). The strain for this 
case may be expressed as 

e^,(r,z) = - ^ /I aJ(1-2v)F^(r.|z|) + | z | ( r , | z | ) ] . 

(-«><z<“), (18) 

where F^ and G^ are also given by (17). 

*The details may be found in [11]. 
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By substituting now from (2), (3)» (4), (16), and (18) into 
(1) we obtain a singular integral equation for the unknown func- 
tion Z(z), Z again has integrable singularities at z=+c,the 
index of the integral equation is +1, and the general solution 
contains one arbitrary constant which is accounted for by the 
si ngl e-val uedness condition of the displacement u^. In the numer- 
ical solution of the problem if one assumes that Z(z) = -Z(-z), 
because of the symmetry properties of the kernels, the single- 
valuedness condition is automatically satisfied [1]. Note that 
in the derivation given in this section the elastic moduli appear 
only through E/E^ . 

3, THE FILAMENT STRESS 


Once the body force Z(z) is obtained the shear stress on the 
filament-matrix interface and the axial stress in the actual fila- 
ment may be obtained by using the general results given in [1]. 
Thus, from (24) or (48) of [1] the interface shear may be ex- 
pressed as 




E^ - E 


Z(z) 


(- c < z < c ) 


(19) 


Similarly the filament stress evaluated by using 

(27) or (49) of [1], which may be expressed as 

c 

/ z(t)dt + Oq + cr^2(b,z) + a^^(b,z) 


* , (-C < z < c) , (20) 

where on the right-hand side of (20) the first term is the stress 
in the auxiliary filament, the second and the subsequent terms 
correspond to the matrix stresses obtained from the problems 

- 10 - 



describBd by Figures 2a -2d. The matrix stresses may either be 
averaged over the area (0<^6 j 0 ^ < 2tt) or be approximated by 

the stress along the filament axis, (r = b, -c<z<c). Follow- 
ing [1] these stresses may be expressed as 

o 

r^ c 


zz 


(b z) = .T° -r / (1 - 2v+ — 2 -)dt 

^ ’ ^ 4 [^2 +(t-z)2]3/2 r2+(t-z)^ 


- 2v + 1^ - ' -^ — J-)] , 

d? + [t-zf 
1 


( 21 ) 


(where d. is given by (5)); 

a^^(b.z) = - a^/7h [Fj,(b.|z|) + | z | G^(b , | z | ) ] , 


( 22 ) 


a^,(b,z) = - /27? I a b FF fb,|z|) + |z|G (b,|z|)] , (23) 

ZZ n = 0 

(with (17), (8), (10) and (15)). 


4. THE STRESS INTENSITY FACTOR 

Because of symmetry on the z = 0 plane “ ^9z " ^ 
matrix, and the stress state in the close neighborhood of and in 
the plane normal to the leading edge of the crack is plane strain. 
Thus the stress intensity factor may be defined by 

K(e) = 11m /27 t( r-a > (r>a, O£0<2ir), (24) 

r^a 

or 

K(e) = - lim /27T(a-r) ^ U2(r,0,+O) , 

r-^a 2(l-v ) 

(0£r<a, 0£e<2Tr), (25) 
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where u^(r,0,+O) is the crack surface displacement on z> 0 side 
of the medium and it is assumed that a = l. From (13. a) and the 
solution given in Section 2 the crack surface displacement in the 
matrix may be expressed as 

U2(r,0,+O) = u^C^j+O) + u^(r,0,+O) 


= ^ /2h [a H fr,0) 


0 0 


+ I a b cos(nN0)H (r,0)] 
n=0 " " " 


(0<r<l, O<0< 2tt) , ( 26 ) 


where 




-az 


Following [10] 


H„(r,+0) = 


^nN,, ^2,1/2 

Jl r(3/2) > (0<r<l), 


(28) 


0 . 


(1 < r < oo) , 

which gives 

2 ^ 

u^(r,e, + 0) = ^ a^b|^r"'^cos(nN0)] , 

(29) 


n = 0 

(0<r<l, 0<6<27 t), 


and 


K(e) 


= _?. [cfQ + I a^b^cos(nNe)] , (0l9< 2ir) . 


(30) 


/F 


n = 0 


5. NUMERICAL RESULTS 

In deriving the integral equation for Z(z) it was pointed 

f 

out in Section 2 (of this paper) that the strain in the auxil- 
iary filament has to be matched with the strain in the matrix 


- 12 - 



along an appropriate generator of the cylinder (<5 = r^, 0£<|i<2Tr, 
-c<z<c), i.e., along the line (<S=rQ, ^“^o’ “C<z<c), The 
effect of this selection on the stress intensity factor for a 
single filament is shown in Figure 3, The three curves given in 
the figure correspond to selecting this generator at the outer- 
most (<J>Q = 0 r*=b+rQ)» at the innermost r’ = b-r^), 

and at an intermediate point (<|)q-tt/ 2 or r’ = b) of the filament 
cross-section. The differences appear to be rather insignificant. 
Hence, for the remainder of the numerical analysis in this paper 
the strains are matched along the line (r = b, -c<z<c). The 
solution of the integral equation is of the form [12] 

Z(z) = f ( z ) / c^ - z^ , (-c<z<c), (31) 

where f(z) = -f(-z) is a bounded function which is obtained 
numerically [13,14]. In the examples given in this paper the 
crack radius was assumed to be the length unit, i.e., a=l. For 
a^l the stress intensity factor ratio given in the figures should 
be read as Kj/ ( 2a^/a/7r ). 

The effect of the modulus ratio E^/E and filament distance b 
in the case of single filament is shown in Figure 4. As expected 
the stress intensity factor decreases for increasing E^/E, de- 
creasing b, and decreasing |e|, where the angle 9 along the crack 
front is measured from the radius going through the filament 
center . 

Figure 5 shows the effect of the filament radius r^ and fila- 
ment length 2c on the angular variation of the stress intensity 
factor for a single filament. 
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Figure 3. The effect of the location of the strain matching 
line on the stress intensity factor. 



Figure 4. The effect of the f i 1 ament-to-matri x modulus ratio 
Ef/E and the f i 1 amen t-to-crack center distance b on the angular 
distribution of the stress intensity factor. 






The effect of the number of filaments, N, is shown in Figure 
6. It may be observed that for N ^4 in the case of b/a = 2 and 
for N^6 in the case of b/a =1.5 (and for c/a = 5, r^/a=0.1, and 
E^/E=100) the stress intensity factor will be independent of 6 
and the problem may be simplified considerably by assuming axial 
symmetry (and by treating the filaments as "shells” of equivalent 
stiffness which are concentric with the crack). 

The effect of the modulus ratio E^/E together with N is 
shown in Figure 7. 

Figures 8 and 9 show the "stress concentration" effect in 

the filament resulting from the existence of the crack. Here, 

is the axial stress in ‘the filament at z = 0 obtained from 

(20) in the presence of the crack for a single filament whereas 

a^^^(O) is the filament stress at z = 0 calculated from (20) in 

b d 

the absence of the crack (i.e., obtained by ignoring and a^^) 
Figure 8 shows the effect of the modulus ratio and the distance b 
Figure 9 shows the effect of the filament dimensions r^ and c, 
and the distance b. 

Figures 10 and 11 show the effect of the modulus ratio E^/E 
on the actual filament stress at z = 0 for a single filament in 
the cracked matrix. The result for E^/E = 1 (i.e., the case of 
no filament) is obtained from the simple penny-shaped crack solu- 
tion (e.g., [4]) as the stress at (r=b, z=0). Note that as 
b/a increases, the filament stress at z = 0 approaches E^/E asymp- 
totically. This is the known result in the uncracked matrix for 
the case of reinforcement by fibers or long filaments. 
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Figure 6. The effect of the number of filaments N on the stress 
intensity factor. 





Figure 7, The effect of the modulus ratio 
intensity factor. 


Ef/Efj on the stress 


1^6 









Figure 9. The ratio of the maximum filament stress in a cracked 
matrix to that in an uncracked matrix for a single filament, the 
effect of filament dimensions. 



I4h 
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Figure 10, Maximum filament stress for a single filament in the 
cracked matrix. 






Figure n. Maximum filame 
cracked matrix, the effect 






Finally, Figure 12 shows the effect of the number of fila- 
ments N and the modulus ratio E^/E on the filament stress at z= 0 
in the presence and in the absence [1] of the crack. The "stress 
ratio" given in the figure is the ratio of the filament stress 

to the filament stress at z = 0 obtained from the solution 
of a single filament in the uncracked matrix. 

The analysis given in this paper may easily be extended to 
filaments with different dimensions, locations, and elastic moduli 
provided they are oriented perpendicular to the plane of the crack 
and this plane is a plane of symmetry with respect to the external 
loads and the geometry of the composite medium. In this case, 
the formulation of the problem would lead to a system of singular 
integral equations in the body forces Z-j (z) , . . . ,Zj^(z) imbedded in 
the filament-matrix interface, N being the number of filaments. 

In the case of symmetrically located large number of filaments, 
the problem may be treated as being axisymmetric with unknown 
body forces imbedded along cylindrical surfaces concentric with 
the crack. These extensions will be considered in the next paper 
which deals with a somewhat more practical problem of infinite 
fibers rather than finite filaments. 
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